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ABSTRACT 



Context. Infrared-dark clouds (IRDCs) are the precursors to massive stars and stellar clusters. GOll. 11-0.12 is a well-studied fila- 
mentary IRDC, though, to date, the absence of far-infrared data with sufficient spatial resolution has limited the understanding of the 
structure and star-formation activity. 

Aims. We use Herschel to study the embedded population of young pre- and protostellar cores in this IRDC. 

Methods. We examine the cloud structure, which appears in absorption at short wavelength and in emission at longer wavelength. We 
derive the properties of the massive cores from the spectral energy distributions of bright far-infrared point sources detected with the 
PACS instrument aboard Herschel. 

Results. We report on the detection and characterization of pre- and protostellar cores in a massive filamentary infrared-dark cloud 
GOll. 11-0.12 using PACS. We characterize 18 cores directly associated with the filament, two of which have masses over 50 Mq, 
making them the best candidates to become massive stars in GOll. 11-0.12. These cores are likely at various stages of protostar 
formation, showing elevated temperature (<T>~ 22 K) with respect to the ambient gas reservoir. The core masses (<M>~ 24 Mq) 
are small compared to that in the cold filament. The mean core separation is 0.9 pc, well in excess of the Jeans length in the filament. 
Conclusions. We confirm that star formation in IRDCs is underway and diverse, and IRDCs have the capability of forming massive 
stars and clusters. 
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1. Introduction 

Star formation occurs in cold, dense molecular cores. In local 
low-mass star-forming regions, this process takes place in rela- 
tively isolated environments and is better understood compared 
to the distant regions in which massive stars form. Specifically, 
it is the identification of the cold precursors to stars - the so- 
called "pre-stellar cores" that has led to a greater understanding 
of the initial conditions of low-mass star formation; however, the 
high-mass counterpart has been more difficult to identify. Thus 
it is not clearly understood how the early stages of massive star 
formation compare to that of low-mass stars. 

Infrared-dark clouds (IRDCs) are complexes of cold (T 

< 20 K), dense (« > 10"* cm"^) molecular gas, now be- 

\ lieve d to be the precursor s to massive stars and star clusters 

. (e.g. iRathborne et al]|2006b . although there is diversity among 

' this class of objects. Early observations revealed that the cold 

dust in IRDCs absorbs the bright Galactic background e mis- 

sion at 8 fim (lEgan et al.1 Il998t iHennebelle et alj 1200 lb and 

emits at sub-mm wavelengths (e.g. ICarev et aP l2000l) . Later, 

high-resolution st udies found that IRDCs exhibit sign ificant 

substructure (e.g. Butler & Tan 2009; Ragan et al.1 120091) . and 

(sub-)millimeter st udies show that they host massive embedded 

star formation (e.g. Rathborne et al.ll2007l; iBeuther & Steinacken 



120071; iHennemann et al.ll2009l) . 

* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. It is open for proposals for observing 
time from the worldwide astronomical community. 



The spectral energy distribution (SED) of massive embedded 
protostars is expected to peak in the far-infrared but is usually 
very faint in the mid-infrared. The coldest cores are not read- 
ily detectable with Spitzer because of sensitivity limitations at 
24 yum and confusion at longer wavelengths due to limited spatial 
resolution. With Herschel (Pilbratt et al. 2010), we now have the 
spatial resolution and the capability to perform large-area maps 
to search for these deeply embedded objects and to characterize 
the SEDs. 

The large, filamentary IRDC GO 1 1.1 1-0.12 was selected as 
a target in the Early Phases of Star Formation (EPoS) Herschel 
key program (PI: O. Kraus e). Together with com plementary 
Spitzer d ata, we use PACS (iPo glitsch et"ai]|2010h and SPIRE 
(iGriffin e t al. 2010) observations to constrain the properties of 
the massive embedded cores within this IRDC. 

Due to its high absorbing contrast against the Galactic back- 
ground emission (particularly at 8 pm), the IRDC G Ol 1.1 1-0.12 
i d — 3 .6 kpc) has been subject to detailed studies. ICarev et all 
1998h confirm the presence of dense molecular gas based on 
H2CO observations. ICarev et aP (l2000l) report 850 and 450 pva 
continuum observations of the thermal dust emission which 
precisely follows the structure seen in absorption in the mid- 
infr ared, and the c olumn density reaches up to 10^^ cm"^ (see 
also iParsons et al.|i2009) . The so-called PI position is seen in 
emission in both 8 and 850 pva images (Fig. 1), a likely indica- 
tor of embedded star formation. 

iPillai et al.1 (|2006a) observe GOl 1.1 1-0.12 as part of an 
ammonia survey, including the (1,1) and (2,2) inversion tran- 
sitions with 40" spatial resolution. The masses of four sub- 
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Fig. 1. The IRDC GOll . 11 -0. 12 from Spitzer to SCUBA wavelengths, (a) Spitzer IRAC 8yum image, (b) Spitzer MIPS 24 pm image. 
The peak SCUBA positions (PI, P2,...) as assigned by Johnstone et al. (2003) are indicated, (c) Herschel PACS 70 pva image, (d) 
Herschel PACS 100 pva image with the positions of analyzed sources (see Sect. 4) indicated. Sources with no 24 pm counterpart 
are circled in cyan, sources with 24 pva counterparts are circled in green, (e) Herschel PACS 160 pva image, (f) Herschel SPIRE 
250 pva image, (g) Herschel SPIRE 350 pva image, (h) 8 pva extinction map (derived from Spitzer data) with (blue) SCUBA 850 pva 
continuum contours (levels are 0.4 0.8 1 .2 Jy beam"'). Source labeling same as in panel (d). The approximate beam size is indicated 
in the lower-left corner of each panel. 



regions (PI, P2, P6, and P7 in the lJohnstone et all(l2003h nomen- 
clature) were found to be several hundred Mq each, though 
other mass-tracing t echniques find highe r values for the total 
filament mass (e.g. [Marshall et al.l l2009h . The NH3 observa- 
tions confirm that temperatures are very low throughout the fil- 
ament, averaging 12 K. In follow-up studies, observations of 
masers dPillai et al v '2006b'). extended excess 4.5 /im emission 
(a tracer of shocked H2 Cyganowski et al. 2008), and molecu- 
lar line emission JPillai et al.ll2006bt IChen et alj|2010bh have all 
supported the presence of ongoing star formation within the fil- 
ament. 



2. Herschel observations and data analysis 

The IRDC c omplex GOl 1.1 1 -0.12 was observed with the PACS 
instrument (Poglitsch et al.' '201(f) on board of the Herschel 
Space Observatory (Pilbratt et al. 2010) on 2009 October 9, 
within the science demonstration program. Scan maps in two or- 
thogonal directions were obtained with the medium scan speed 
of 20"/s. Customized scan leg lengths (600", 1200") and scan 
position angles (12.5°, 102.5°) were employed in order to adapt 
for the elongated shape and orientation of the science object. In 
total, 2280, 2280, and 4560 seconds were spent for scanning on 
target in the blue (70 pm), green (100 pm), and red (160 pm) 



PACS filters, re spectively. The raw data were reduced with the 
HIPE software (IOttll2010l) . version 3.0, build 455. Besides the 
standard steps leading to level- 1 calibrated data, a second-level 
de-glitching as well as a coiTection for offsets in the detector 
sub-matrices were performed. Finally, the data were highpass- 
filtered using a large median window of 600", in order to remove 
the effects of bolometer sensitivity drifts and 1/f noise along the 
course of the data acquisition. 

Since we are mainly interested in the population of un- 
resolved point sour ces, we used the IDL Starfinder program 
jDiolaiti et al.l l200u !) to identify point sources and apply PSF 
photometry to the PACS data. Reference PSFs (taken on the as- 
teroid Vesta) were provided for all 3 filters by the PACS instru- 
ment team. 

Maps at 250, 350, and 500 pm were obtained with SPIRE 
(Griffin etal. 2010) on 19 October 2009. Two 720" scan legs 
were used to cover the source. Four repetitions resulted in 340 s 
of scanning time with the nominal speed of 30"/s. The data were 
processed within HIPE (Ott 2010) with the standard photometer 
script. Since for these observations no cross-scan data were se- 
cured, the resulting maps still show residual stripes along the 
on e scan direct i on. T he iterative de-striping algorithm devised 
bv lBendo et alJ (I2OIOI) was invoked to mitigate this effect. 
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Table 1. Photometric Properties of Protostars 



Src. 




RA 


Dec. 


Temp. 


Lum. 


Mass 


Notes 


Num. 


(J2000) 


(J2000) 


(K) 


(Lo) 


(Mo) 




1 


8 


10:47.0 


-19:10:09 


21 


6 


3 




2 


8 


10:40.3 


-19:11:34 


23 


26 


9 


P3,e 


3 


8 


10:39.9 


-19:16:39 


18 


5 


5 




4 


8 


10:38.8 


-19:19:01 


24 


8 


1 


e 


5 


8 


10:38.2 


-19:19:50 


25 


66 


10 


e 


6 


8 


10:45.9 


-19:20:46 


16 


15 


40 




7 


8 


10:46.3 


-19:21:46 


20 


6 


4 




8 


8 


10:33.7 


-19:22:08 


21 


25 


9 


P2, e 


9 


8 


10:28.4 


-19:22:29 


24 


1346 


240 


Pl,e 


10 


8 


10:03.8 


-19:23:06 


23 


27 


7 


e 


11 


8 


10:16.4 


-19:24:19 


19 


9 


7 




12 


8 


10:19.4 


-19:24:24 


21 


11 


5 


e 


13 


8 


10:09.2 


-19:25:55 


26 


23 


2 


P2,e 


14 


8 


10:06.9 


-19:26:00 


20 


9 


5 




15 


8 


10:05.5 


-19:26:34 


24 


58 


11 


P7, e 


16 


8 


10:04.7 


-19:27:43 


24 


25 


5 


P6, e 


17 


8 


10:04.0 


-19:28:16 


23 


28 


6 


P6 


18 


8 


10:07.4 


-19:29:15 


23 


139 


82 


P6,e 


19 


8 


10:25.1 


-19:17:49 


24 


478 


79 


e 


20 


8 


10:27.1 


-19:18:12 


23 


45 


11 


e 


21 


8 


10:07.7 


-19:18:39 


21 


19 


9 




22 


8 


10:42.7 


-19:19:58 


21 


11 


5 


e 


23 


8 


10:03.8 


-19:21:04 


20 


8 


5 


e 


24 


8 


10:13.8 


-19:22:04 


17 


9 


15 


e 



N otes - "P" lab els indicate spatial correspondence to regions identified 
by lJohnstone e t al. (2003). "e" indicates where we infer the presence of 
an embedded protostar based on a 24 /jm detection. Objects above the 
horizontal line are located "on" the IRDC filament; sources below are 
"off." 



3. Results and discussion 

In Spitzer images, GOll.11-0.12 appears mostly dark at both 
8 and 24 fim, as shown in Fig. [T] though there are some in- 
dications of point sources at 24 /im. However, the nature of 
these objects remains undefined because they lack counterparts 
at shorter wavelengths. The IRDC filament is dark in the PACS 
70 and 100 fim filter images, but shows little evidence of absorp- 
tion at 160 i-im. In the continuum maps, st arting with 250 ^m 
SPIR E map to the SCUBA 850 //m map dPi Francesco et alJ 
12008 ). shown in contours in Fig.[Tl the extinction filament seen 
at shorter wavelengths now appears in emission. 



3.1. Characterization of embedded cores 

Herschel provides both the ability to make large-scale maps and 
the spatial resolution to isolate point sources from bright back- 
ground emission for the first time. At the given scales between 
18000 and 43000 AU (our resolution limit from the 5 to 12" 
PACS beam), we probe at the size scale of "cores," as defined by 
iBergin & Tafallal (J2007). We detect cores seen as weak 24 /im 
point sources and new cores with no 24 ircn. counterparts (see 
Fig. 1). The point sources follow the general structure of the 
extinction/emission maps quite well. However, there are a few 
objects associated with bright emission to the north and west of 
the extinction map which are potentially associated with a more 
evolved stellar cluster. 

The PACS observations cover a crucial portion of the spec- 
tral regime in which cold cores are luminous. For point sources 
with detections in the PACs bands (70, 100, and 160 //m), we 
fit the SED with modified Planck blackbody functions account- 
ing for the wavelength dependence of optical depth assuming 



iDraine & Led (Il984t) dust composition. Note that each wave- 
length has a different beam size, thus the flux included in our 
photometric measurements is taken over correspondingly dif- 
ferent spatial scales. We do not include the 24 //m data point 
because the emission becomes optically thick, negating the as- 
sumptions of the optically thin blackbody fit. The SPIRE data 
are important to constrain the large-scale properties of the fila- 
ment, however with the poorer spatial resolution, we do not in- 
clude them in the blackbody fits because the PACS cores cannot 
be identified. The SCUBA data trace the colder surrounding fil- 
ament, for which the peaks are not necessarily coincident with 
PACS point sources, so these data are also excluded from the 
SED fit. 

In Table [T] we present the position, model temperature, lu- 
minosity, and mass from the single-component blackbody model 
of the point sources with detections in all three PACS bands. 
The mean mass, luminosity, and temperatures of the sources are 
24 Mq, 96 Lq, and 22 K, respectively. Assuming fiducial uncer- 
tainties for the PACS band fluxes (10%, 20%, and 20% for 70, 
100, and 160 yum respectively), the temperature, luminosity, and 
mass values presented are accurate to approximately 4%, 15%, 
and 40%, respectively. We also note in Table [1] which sources 
have 24jum detections (16/24). 

Cores with a 24 jum detection likely host a protostar. Cores 
with no 24 /vm detection may be prestellar, but geometry may 
limit detectability at 24 /vm (e.g. obscuration by dusty disk). 
There are an additional 20 point sources not reported here de- 
tected in only two of the PACS bands that are also candidate 
prestellar cores. Sources with no 24 jum detection (and even 
those with no 70 ixn\ counterpart) ca nnot be definitively clas- 
sified as prestellar. IChen et alJ (1201 Oa) has shown that (in a low- 
mass core) the absence of 70 /im emission cannot rule out the 
presence of a protostar. Further observations (e.g. outflow signa- 
tures) are needed to verify the absence of deeply embedded pro- 
tostars. Perhaps the best candidates for truly "prestellar" objects 
are the SCUBA peaks without any PACS emission. These are 
column density peaks at an even more primitive stage, however 
without interferometric follow up at comparable resolution to the 
Herschel observations, the substructure of the prestellar material 
can not be fu lly characterized. For exa mple, iHenneman n et all 
(.2009.) and .Beuflier & Hennind (l2009h have shown that such 
dark regions in IRDCs fragment into a few massive cores. 

In all, 18 of the 24 cores (75%) are located "on" the IRDC, 
which we define as positions detected in the SCUBA map or on 
the 8 jum extinction map (see Fig. 1) in regions where SCUBA 
data do not exist. There are two protostellar cores on the filament 
which have masses larger than 50 Mq, making them excellent 
candidates to host massive protostars. We detect the PI source 
(9 in Table 1) at all wavelengths, and we estimate the core mass 
to be 240 Mq, the highest in our sample. 

The large scale map provides much information about indi- 
vidual parts of the filament, but now with high-resolution maps, 
we explore some detailed properties of subregions. As an exam- 
ple, we discuss a region located south and west of PI, which 
includes sources 11 and 12. In Fig.|2l we show the 24 /vm and 
70yum image of this small region spanning 1.7 pc, and the SEDs 
of the two sources. The masses of cores 11 and 12 are 7 and 
5 Mq, respectively, and the total mass reservoir estimated in the 
region is roughly 700 Mq (derived from SCUBA 850 //m dust 
continuum). The western source (11, bottom plot) is not detected 
at 24 jum (a pre-/protostellar core), and the eastern source (12, 
top plot) has a 24 /im detection in excess of the SED fit to the 
PACS data (a protostellar core). While both in early evolution- 
ary stages, we infer that source 12 is more evolved than source 
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Fig. 2. Left: SCUBA contours (levels: 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 
0.7, 0.8, 0.9 Jy beam"') plotted over MIPS 24 /jm image (top) 
and PACS 70 fim image (bottom). Arrows indicate the location 
of point sources 11 (east) and 12 (west) from Table 1. Right: 
Spectral energy distributions of sources 1 1 (bottom) and 12 (top) 
detected within the extincted region. 



11, even though it is likely that they originate from the same 
natal clump. We also find that the local SCUBA emission peak 
does not directly coincide with either point source (see Fig. |2]i, 
with the nearest source (11) offset by 19" from the peak. The 
core masses are small compared to the local clump mass within 
the filament. Interferometric follow-up is key in determining the 
substructure within th e region and the full cen sus of pre- and 
protostellar objects (cf. lHennemann et al.ll2009t) . 

The SEDs presented in Fig. |2] include the PACS data used 
in the blackbody fit, limits on the Spitzer flux at 3.6, 4.5, 5.8, 
8.0, and 24 yum, except in the case of core 12 where there is a 
24 //m detection. The data point at 850 yum on the SEP is the 
flux estimated within a 23" SCUBA beam dPi Francesco et al] 
[2008) at the position of the point source (corresponding to 83000 
AU). In all cases where there are SCUBA data, the flux within 
the SCUBA beam is in excess of the SEP predicted from the 
blackbody fit to the PACS data. This is not surprising because of 
the larger aperture and that the SCUBA measurement probes the 
emission from the colder filament, which includes only a small 
contribution from the war mer (~23 K) core. Assuming a fila- 
ment temperature of 1 2 K (iPillai et al.lf2006ah and k^^q - 0.01 
cm^g ' (Praine & Lee 19841) . the flux within the SCUBA beam 
corresponds to a mass 40 - 80 times the core mass. The SCUBA 
beam probes the mass in the outer envelope of the cores. We do 
not claim to constrain this envelope size, but rather provide a 
uniform benchmark against which we can compare the cores to 
their immediate surroundings. 



3.2. The nature of star formation in IRDCs 

The small aperture of Herschel allows us to probe the embed- 
ded population on the "core" scale in GOl 1.1 1-0.12. The mean 
core separation is 0.9 pc, which is many times the Jeans length, 
Ry (for T = 12 K and 10"* - 10^ cm-^ R,, = 0.06 to 0.2 pc), in- 
dicating that Jeans fragmentation seems not to be dominant in 
this IRPC. Cores with a 24 //m detection are likely protostellar 
cores, and cores with only PACS detections can be either prestel- 
lar or protostellar ("pre-/protostellar"). The SEP fits give core 
masses from 1 to 240 Mq, showing that the core population, and 
therefore the star formation activity in this IRPC, is diverse. The 
core temperatures are fairly homogeneous throughout the fila- 
ment, averaging 22 K, which is ~10 K warmer than the gas tem- 



perature measured in the large-scale ammonia maps. This may 
indicate the presence of embedded star formation activity that 
would not be dete ctable with low-resolution observations (e.g. 
iPillai et al.ll2006ah . 

The cores presented in this study represent a key early stage 
in protostar formation. This important phase has until now been 
inaccessible, since Herschel allows large, high-resolution sur- 
veys of the far-infrared wavelength regime for the first time. 
There are certainly also more primitive phases of star forma- 
tion within IRPCs. The filament is comprised of thousands of 
Mq of cold gas, most of which is not in cores detectable with 
Herschel PACS. The abundant cold gas may serve as a reservoir 
from which the cores accrete mass. 



4. Conclusion 

Herschel enables us to perform large-scale, multi-wavelength 
surveys of star formation regions with unprecedented sensitiv- 
ity and efficiency in the far-infrared. In this contribution, we re- 
port the detection of embedded pre-/protostellar cores in the fil- 
amentary IRPC GOl 1.1 1-0.12. We characterize the SEPs of 24 
cores. The average mass, luminosity, and temperature of all the 
sources is 24 Mq, 96 Lq, and 22 K. Compared to the mass of 
the filament estimated from other surveys, there is an immense 
reservoir of mass that is not confined to pre-/protostellar objects. 
The cores have an elevated temperature compared to that esti- 
mated from the IPillai et al.l (2006a) large scale ammonia survey 
of the filament, which may be due to localized heating of the gas 
by protostars. 

In forthcoming work, we intend to extend this investiga- 
tion to perform more detailed radiative transfer modeling of the 
SEPs, connect the forming protostars to the kinematics, and ex- 
tend these analyses to a larger sample of IRPCs observed as part 
of the EPoS program. 
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